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Abstract. - We demonstrate the existence of superscarring in the LiNC^LiCN isomerization 
reaction described by a realistic potential interaction in the range of readily attainable experimen- 
tal energies. This phenomenon arises as the effect of two periodic orbits appearing "out of the 
blue" in a saddle-node bifurcation taking place in the dynamics of the system. Potential practical 
consequences of this superlocalization in the corresponding wave functions are also considered. 
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Introduction. — The localization of certain quantum 
eigenstates along unstable periodic orbits (PO) in classi- 
cally chaotic systems - i.e., quantum scarring [1] - has 
profound consequences for the dynamics of quantum sys- 
tems. It was first observed in numerical calculations on 
the Bunimovitch stadium billiard [2]. At the most ba- 
sic level, scarring [3, 4] is a spectacular example of the 
influence of classical structures, i.e. PO [5] and their as- 
sociated manifolds [6-8] on quantum behavior [9]. Scar- 
ring is responsible, for example, for significant deviations 
from the predictions of random matrix theory [10]. It has 
also been conjectured recently that they are important in 
open systems to understand the semiclassical structure of 
resonances [11, 12], the number of which has been seen 
to scale with h according to a fractal Weyl law [13, 14]. 
At the same time, scarring has important practical con- 
sequences in various fields, including tunnel currents in 
mesoscopic devices [15], quantum dots [16], optical cavi- 
ties [17], optical fibers [18], and very recently the existence 
of relativistic quantum scars has been demonstrated [19] 
in the emerging field of graphene applications [20] . 

The vast majority of the work published on scars has 
concentrated on fully chaotic systems, where the POs 
are isolated and unstable. However, Keating and Prado 
(KP) [21] recently considered the specific peculiarities of 
generic systems (with mixed dynamics), where bifurca- 



tions of the scarring PO occur. By extending Bogomolny 
theory [3] , they found that the density localization around 
bifurcated POs are wider and with amplitudes at least as 
large as those far from bifurcations. This led KP to in- 
troduce the term superscars to name this new localization 
phenomenon which they illustrated with numerical exam- 
ples in perturbed cat maps. The influence of bifurcating 
orbits on the density of states was studied by Schomerus 
and Sieber [22] , who also found different scaling exponents 
on the dependence of such magnitude with h. Due to the 
nature of the arguments used by KP, which only apply 
to high excitation energies (K — ► limit), an open ques- 
tion has remained and concerns the transferability of KPs 
conclusions beyond maps to actual physical systems over 
experimentally realizable energy ranges. 

In this Letter, we address this crucial point, namely 
the existence of a universal behavior in the exponents of 
the moment distribution characterizing superscars [21]. 
This is done in the context of a realistic model for the 
LiNC^LiCN isomerization reaction [23]. More specifi- 
cally, for this system we investigate the superscarring as- 
sociated with POs which originate in a saddle-node (SN) 
bifurcation of the isomerization dynamics [24] . 

Isomerization plays a fundamental role in many bio- 
chemical processes and it is an ideal benchmark for dy- 
namical systems theory [23,24]. Also, the example that 
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we have chosen can be considered as a prototype for this 
kind of processes (see description of the model below), 
which has the added bonus of exhibiting a rich and com- 
plex vibrational dynamics. Moreover, the superscarring 
phenomenon discussed here is closely related to some par- 
ticular aspects of the transition state (TS) theory (TST) 
[25] at quantum level. This theory, that was developed 
in the 1930's after the seminal work of Eyring [26] and 
Wigner [27], envisages the physical processes associated to 
chemical reactions as "skiing the reaction slopes" [28] or 
landscapes [29] , so that the barriers existing between reac- 
tants and products are surpassed. Quantum mechanically, 
these barriers accommodate the corresponding Gamow- 
Siegert resonant states [30,31], and the topology of the 
associated wave functions configures the doorways con- 
trolling the quantum reactivity [32] . 

In the title reaction there exists, apart from the reaction 
barrier corresponding to the saddle point in the potential 
energy surface that was previously studied in Ref . [33] , an- 
other interesting barrier of dynamical origin [34,35], asso- 
ciated to the POs arising from a SN bifurcation in the cor- 
responding dynamics [24]. According to KP's results the 
characteristics of this barrier changes dramatically with 
the excitation energy. At high energy, where the two POs 
are isolated the quantum density localization effect due 
to scarring is less pronounced. On the contrary, when 
lower values of the energy, close to the bifurcation point, 
are considered a strong superlocalization takes place. The 
fundamental difference in the widths taking place in this 
two different situations is expected to have some relevant 
influence in the corresponding flux [36], leading to a dif- 
ferent behavior in the respective reactive dynamics. 

Finally, it should be pointed out that other interesting, 
although completely different, effects of SN bifurcations 
in the vibrational dynamics of molecular systems (HCP 
molecule) have been described in the literature [37]. 

Model and calculations. — As a working model, 
we choose the LiNC/LiCN isomerizing molecular system, 
which has been extensively studied in connection with 
quantum chaos [23]. This system is representative of a 
large class of small polyatomic molecules, which exhibit 
similar behavior, mainly due to the presence of a large am- 
plitude (floppy) motion in one of the vibrational modes. 
This class includes other cyanides, such as HCN/HNC 
[38]), alkaline cyanides [39], and other similar especies, 
such as, HCP [37], the HO2 radical, or the van der Waals 
complexes. 

Our system can be modelled by a realistic two degrees of 
freedom vibrational (J = 0) Hamiltonian, which is given 

by 

Here we use Jacobi coordinates, which are common for 
example in scattering theory [40] and celestial mechan- 
ics [41]. In our case these coordinates, (r, R, ■&), are de- 




Fig. 1: Definition of the Jacobi coordinates used in our calcu- 
lations. 

fined as the C-N distance, the distance between the Li 
atom and the center of mass of the CN fragment, and 
the angle formed by these two vectors, respectively (see 
Fig. 1). The corresponding reduced masses are defined as 
Hi = m Li m CN /m LiCN and fi 2 = m c m N /m CN . The C-N 
motion is kept frozen at its equilibrium value, r e , since 
the associated frequency is high, thus decoupling very ef- 
fectively from the rest of coordinates in the molecule. The 
interaction potential is expressed as a series of Legendre 
polynomials, 

9 

V(R,d) = '£ i Px(coa#)vx(R), (2) 

A=0 

whose coefficients consist of a sum of short-range and long- 
range contributions whose explicit forms have been taken 
from the literature [42]. It is shown in the top panel of 
Fig. 2 as a contours plot. As can be seen, it presents two 
wells at i3 = and 180°, corresponding to the linear iso- 
mers LiCN and LiNC, which are separated by a modest 
energy barrier. The motion in 1? is very floppy and chaos 
sets in at moderate values of the excitation energy. Su- 
perimposed to the potential surface, we also present in 
the plot the POs which are relevant for our discussion on 
superscars, at two different values of the interaction en- 
ergy. They all run almost vertically, and correspond to 
a 1:1 resonance which appear "out of the blue" due to 
a SN or tangent bifurcation at £bif = 3440.64 cm -1 , as 
discussed in Ref. [24]. At this energy, the two orbits are 
indistinguishable, since they both coalesce in the central 
(shorter) trajectory shown in the figure. As energy in- 
creases the two POs separate, the unstable one moving to 
the right and the stable towards the left. This stability 
is later lost for E > 4960.1 cm -1 . These two POs have 
also been plotted in the top panel of Fig. 2 for an energy 
far away from the bifurcation point, namely E\ = 9196.0 
cm -i ~ 3.Ebif- 

The corresponding classical dynamics can be adequately 
followed by computing Poincare surfaces of section (SOS) 
using the minimum energy path, i? e (i?), connecting the 
two isomers, as the sectioning plane. This choice requires 
an extra canonical transformation, 
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Fig. 2: (Top) Potential energy surface for LiNC/LiCN. The 
minimum energy path (dotted line), and three saddle-node 
periodic orbits (full line) at _Ebif=3440.64 cm -1 (central orbit) 
and _Ei=9196.0 cm -1 (lateral trajectories) relevant to our work 
have been plotted superimposed. (Middle and bottom) Com- 
posite Poincare surface of sections for the same energies. The 
position of the fixed points corresponding to the periodic orbits 
are represented with full squares. 

to make the SOS an area preserving map [43]. Composite 
SOSs for Ebif and E\ are shown in the middle and bottom 
panels of Fig. 2, respectively. Two comments are in order. 
First, the dynamics is mixed, being chaotic in the region 
of interest - i.e. where the two bifurcated POs sit - and 
regular around the isomer configurations. Second, the ob- 
vious signature of a cantori [44] , seen as an accumulation 
of points next to the region of stability around $=180°, is 
apparent for £?bif> and this bottleneck disappears at E\. 

As mentioned before, one critical element in our work 
is the computation of the scarred wave functions. The 
reason for this should be clear. Since we want to perform 
our study far from the high E (low K) semiclassical limit, 
where the quantum oscillations in the local configurations 
and energy averages implied in scarring [3] easily cancel 
out, an algorithm able to introduce the subtle dynamical 
effects associated to the scarring PO in the scar functions 
is forcefully needed. In our case, this is achieved by using 



the method described in [45] consisting in time evolving a 
wave packet, \(f>(t)}, initially launched at one of the turning 
points of the PO under study, obtaining afterwards all the 
relevant information by Fourier transform. 

We start with the (infinite resolution) stick spectrum: 

/+oo 
dt e iS */ ft <0(O)|0(t)> 
- CO 

= ^|(0(O)|n)| 2 5(i?-i? n ), (4) 

where kets \n) represent the eigenfunctions of the system. 
An example of the results obtained from such calculations 
is shown in Fig. 3. As can be observed, the distribution 
of sticks, both in separation and length, is rather irregu- 
lar, but when examined closely they are seen to be nicely 
grouped in clumps. This regularity is due to the recur- 
rences in the autocorrelation function induced by the PO 
motion [5], and it shows up as well-defined bands in the 
low resolution version of the spectrum, 

dte -t 2 /T 2 e iEt/H {mm) , (5) 

-oo 

where the cutoff parameter T has been introduced in or- 
der to eliminate the long term dynamics not directly as- 
sociated to the scarring process. The result for T equal 
to the period of the scarring PO is displayed in Fig. 3. 
Here, a series of bands are seen to appear centered at the 
discrete values of the energy, E^° , that are obtained semi- 
classically using the suitable Bohr-Sommerfeld (BS) quan- 
tization condition of the action along the PO [see eq. (8) 
below]. In a second step, we focus our attention in one sin- 
gle of these bands, and calculate the scar wave function, 
\ip n ), associated to it by Fourier transforming the original 
packet, \<f>(t)}, at the corresponding BS energy, E® s , 

/oo 
dte-^/^e^-^j^o)). (6) 
-oo 

Here we choose the Ehrenfest time [46] as the value for the 
cutoff parameter T, since it is known [7] that this proce- 
dure guarantees that the dynamical information concern- 
ing the hyperbolic structure associated to the scarring PO 
is embedded in the topology of \ip n )- By expanding \4>(0)) 
in the basis set of the system eigenstates, the integral in 
cq. (6) can be solved explicitly, and one obtains 

|Vn> = ^T E ^- (E -- EBS)2T2/Ah2 {n\m)\n)- (7) 

n 

This equation clearly shows how the scarred states defined 
by us fulfill the two most important criteria accepted in 
scar theory. On the one hand, they are built by averag- 
ing over the eigenstates within an energy window of width 
K/Te around E® s , in the spirit of Bogomolny's interpre- 
tation of scarring [3]. On the other hand, the particular 
average that we use incorporates the dynamical correla- 
tions existing in ((p(0)\(p(t)) , in agreement with the view 
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Fig. 3: Infinite resolution spectrum and its low resolution ver- 
sion generated from a wave packet localized along the central 
1:1 periodic orbit on the top panel of Fig. 2. The inset shows 
the scar function ip3, whose mean energy roughly correspond 
to E hi{ . 



introduced by Heller [5]. It should be stressed that these 
correlations are so strong that the scar functions acquire 
a very pronounced semiclassical character, independently 
of particular details of the individual eigenfunctions used 
in the construction, and in this sense our method is very 
robust. In the inset to Fig. 3 we show as an example \ipa), 
whose (mean) energy roughly corresponds to E^n and it 
is constructed with the eigenstates in the fourth clump. 
Let us imagine now that h were reduced to one half of its 
value. The energy density would increase, and the clumps 
of Fig. 3 would shift to the left, i.e. to the region of smaller 
energies. The eigenfunctions contributing to each clump 
would then be different, but our computed scar wave func- 
tions would be the same. 

In conclusion, a scar function is only characterized by 
the number of excitations along the PO, and this number 
depends both on the values of E and h, which are given 
by the Bohr-Sommcrfcld quantization condition. Another 
relevant characteristic of the scar functions is related to 
the direction transverse to the PO. In the case of an unsta- 
ble PO, a hyperbolic Hamiltonian dominates the flux in its 
neighborhood and the wave packet, 1 0(0)) , moves quickly 
along the associated manifolds. On the other hand, at 
bifurcations the dynamics acquire some stability, and the 
original wave packet move away from the central PO, but 
in a more slowly way. These two different behaviors have 
an influence in the width of the related scar functions, as 
will be demonstrated in the present article. 

Results and discussions. — Let us next quantify the 
scarring effect on ip n due to the 1:1 SN POs that we are 
considering. As described by KP [21], this can be achieved 
with the aid of the second moment of the wave function 
amplitude distribution averaged over regions shrinking as 
h — > but containing an increasing number of do Broglie 
wavelengths. In the case of an isolated unstable PO this 



magnitude follows, in the limit, a universal power law h 1 ^ 2 . 
However, for just bifurcated POs of the type that we are 
considering here, i.e. tangent bifurcation, this expression 
goes as ft 1 / 3 due to the superscarring effect of the bifur- 
cated orbits [21]. In order to carry out this calculation, we 
first select two values of the vibrational energy, one cor- 
responding to the bifurcation energy, E^f, and the other 
one, Ei, sufficiently far away from it to ensure that the PO 
is isolated. Second, we compute a series of scar functions, 
\ip n ), with different values of the excitation quantum num- 
ber n at the same energy, E. In order to do this, we vary 
H according to the BS quantization condition for the PO 



h{n)=S E (n+ V -f) \ 



(8) 



where 5_E bif =3.5, VE bl[ = 2, Se^ =13.75, ve 1 = 3 in our 
case. Notice that for the second energy the Maslov index 
has increased by one due to the appearance of a self-focal 
point in the PO. The net effect of the H scaling represented 
by (8) is to shift the center of the ra-th band to the value of 
the working energy, E^n or Ei . Some scar wave functions 
computed in this way are shown in Fig. 4. 

As can be seen, they all correspond to different excita- 
tions in the R motion with no nodes in the $ coordinate. 
Moreover, their widths (transverse to the PO) change very 
slowly at Ebif, which is a distinctive property of the phe- 
nomenon of superscarring. In order to make this assertion 
quantitative, we evaluate the variance of the probability 
density at a fixed value of R, 



4(n) = J dtf 



if 



W(R) 



1 



W{R) 



(9) 

with W(R) = J dti \ip n (R,ti)\ 2 . To illustrate this proce- 
dure, the profile of each function in Fig. 4 is shown for 
R = 4.2 a.u. in the insets. Then, we calculate the mean 
width weighted by W(R) as 



a(n) 



I 



dR W(R)a 2 R (n) 



"I 1/2 



(10) 



The corresponding results for E^f and Ei are shown in 
Fig. 5. Several comments are in order. First, as can 
be seen in the log-log plot both sets of data follows re- 
markably well a scaling power law, (n + z^;/4)~ Q , with 
the excitation quantum number, n (equivalent to K). Sec- 
ond, from the slopes of least squares fittings we have that 
a Ebit = 0.32 ± 0.01 and a El = 0.53 ± 0.02. Third, these 
values agree extremely well, within the numerical error, 
with the universal behavior predicted in the theory of KP, 
according to which the exponent is 1/3 at the bifurcation 
point, and 1/2 far from the bifurcation, respectively. 

This last result confirm beyond any doubt that the SN 
POs studied in this Letter, which corresponds to an ex- 
cited molecular motion for the LiNC/LiCN isomerizing 
molecular systems taking place mainly along the stretch- 
ing coordinate, R, is able to produce superscarring on 
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Fig. 4: Some scar functions computed with Eq. (6) at different values of h obtained from Eq. (5). The scarring periodic orbit 
is superimposed. The insets show the wave function profiles along one of the cuts at constant R (shown as short horizontal 
segments) which are used in Eq. (9). 
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Fig. 5: Log-log plot of the averaged scar wave functions 
widths, computed by Eq. (10), vs (n + i> E /4), for two values 
of the vibrational energy: Ebu = 3440.64 cm— 1 (open circles) 
and Ei = 9196.0 cm— 1 (full squares). Also shown are the 
least squares fitted straight lines corresponding to power laws, 
(n + u E /4 : y a with scaling exponents 0.32±0.01 and 0.53±0.02, 
respectively, indicating that we are in the presence of super- 
scarrings. 



Conclusions. — Summarizing, in this Letter we have 
reported for the first time conclusive evidence for the exis- 
tence of superscarring associated to bifurcated POs [21] in 
a realistic molecular system. More interestingly, this hap- 
pens in the regime of very low excitations along the PO, 
which opens the possibility for considering the influence of 
scarring in quantum processes where the probability flux 
is strongly affected by the existence of dynamical barri- 
ers. In particular, the effect of scarring on LiNC^LiCN 
isomerization rates and specially its dependence on the 
value of the vibrational excitation energy is a topic that 
is currently under investigation. 
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the corresponding quantum wave functions for energies 
near the bifurcation point, and this effect disappears as 
we move away form this point. 
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